Phase gate of one superconducting qubit simultaneously controlling n qubits in a 
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We propose how to realize a three-step controlled-phase gate of one superconducting qubit simul- 
taneously controlling n qubits selected from iV qubits in a cavity (1 < n < N). The operation time 
of this gate is independent of the number n of qubits involved in the gate operation. This phase 
gate controlling at once n qubits is insensitive to the initial state of the cavity mode and can be 
used to produce an analogous CNOT gate simultaneously acting on n qubits. 
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Introduction. — Quantum information processing has 
attracted considerable interest during the past decade. 
The building blocks of quantum computing are single- 
qubit and two-qubit logic gates. So far, a large num- 
ber of theoretical proposals for realizing two-qubit gates 
in many physical systems have been proposed. More- 
over, two-qubit controlled- not (CNOT) or controlled- 
phase (CP) gates have been experimentally demonstrated 
in, e.g., cavity QED [![, ion traps 0, NMR quantum 
dots [J| , and superconducting qubits [B|, @] . 

Attention is now shifting to the physical realization of 
multi-qubit controlled gates (e.g., Q) instead of just two- 
qubit gates. It is known that multi-qubit controlled gates 
play a significant role in constructing network quantum 
computation circuits. When using the conventional gate- 
decomposition protocols to construct a multi-qubit con- 
trolled gate [1, the procedure usually becomes com- 
plicated (especially for large n), as the number of single- 
qubit and two-qubit gates required for the gate imple- 
mentation heavily depends on the number n of qubits. 
Therefore, building a multi-qubit controlled gate may be- 
come very difficult since each elementary gate requires 
turning on and off a given Hamiltonian for a certain pe- 
riod of time, and each additional basic gate adds experi- 
mental complications and the possibility of more errors. 

Several methods for constructing phase gates with n- 
control qubits acting on one target qubit based on cavity 
QED or ion traps have been recently proposed hJ 11,12, 
|l3j. These methods open a new way for realizing quan- 
tum controlled-phase gates with multiple control qubits. 
However, we note that these proposals [l(J U, 12, HH 
cannot be extended to perform a different type of sig- 
nificant multi-qubit controlled-phase gate, i.e., quantum 
controlled-phase gates with one qubit controlling n tar- 
get qubits. In this work, we propose how to realize a 
three-step controlled-phase gate of one superconducting 
(SC) qubit simultaneously controlling n qubits selected 
from N qubits in a cavity (1 < n < N). To achieve 
this, we construct an effective Hamiltonian which con- 
tains interaction terms between the control qubit and 
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FIG. 1: (a) A controlled-phase (CP) gate simultaneously act- 
ing on n target qubits (2, 3, n + 1), i.e., n two-qubit CP 
gates. Here, Z represents a controlled-phase flip on each tar- 
get qubit. Namely, if the control qubit (i.e., qubit 1) is in the 
state |1), then the state |1) at each Z is phase-flipped as |1) 
— > — 1 1), while the state |0) remains unchanged, (b) CNOT 
gate simultaneously controlling n qubits (2, 3, n + I), ob- 
tained from our n-qubit CP gate. The symbol © represents 
a CNOT gate on each target qubit. If the control qubit is in 
the state |1), then the state at © is bit flipped as |1) — > |0) 
and |0) — » |1). However, when the control qubit is in the 
state |0), the state at © remains unchanged, (c) Diagram of 
a superconducting (SC) charge qubit. (d)N SC qubits are 
placed in a microwave cavity, from which a subset of qubits, 
selected for the gate, are coupled to each other via the cav- 
ity mode. In (a) and (b), the qubits (n + 2, n + 3, N) 
are not involved in the gate operation, by setting their $ = 
$o/2, Vt = e/C g , and F 9 ac = 0. 



each subordinate or target qubit. We will denote this 
n-target-qubits control-phase gate as a NTCP gate (see 
Fig. 1(a)). As shown below, our proposal has the fol- 
lowing advantages: (i) the n two-qubit controlled-phase 
gates involved in the NTCP gate can be performed simul- 
taneously; (ii) the operation time required for the gate 
implementation is independent of the number n of qubits 
involved in the gate operation; (iii) this approach is in- 
sensitive to the initial state of the cavity mode, and thus 
no preparation for the initial state of the cavity mode is 
needed; (iv) no large detuning between the qubits and 
the cavity mode is required and thus the gate operation 
can be speeded up; and (v) the proposal is remarkably 
simple, requiring only three basic operations. Note that 
a CNOT gate of one qubit simultaneously controlling n 
qubits, shown in Fig. 1(b), can also be achieved using 
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the present proposal. This is because the n-target-qubits 
CNOT gate is equivalent to a NTCP gate plus a single- 
qubit Hadamard gate acting on each target qubit before 
and after the NTCP gate. To the best of our knowledge, 
our proposal is the first to demonstrate that a power- 
ful phase gate, synchronously controlling n qubits, can be 
achieved with superconducting charge qubits in a cav- 
ity, which can be initially in an arbitrary state. This 
proposal is quite general and can be easily extended to 
other physical systems (such as atomic qubits, quan- 
tum dots, and superconducting flux or phase qubits), 
since the relevant effective Hamiltonian [Eq. (fTO]) below] 
can be constructed by applying suitable external driv- 
ing pulses. We believe that this work is of general in- 
terest and significance because it provides a simple pro- 
tocol for performing controlled-phase (or controlled-not) 
gates with multiple-target-qubits, which are important 
in quantum information processing such as entanglement 
preparation [3], error correction [lH, quantum algo- 
rithms (e.g., the Discrete Cosine Transform [l6j]), and 
quantum cloning [l7j |. 

Model. — The superconducting charge qubit consid- 
ered here, shown in Fig. 1(c), consists of a small box, 
connected to a symmetric superconducting quantum in- 
terference device (SQUID) with capacitance Cjo and 
Josephson coupling energy Ejq. In the charge regime 
A > E c > Ej > k B T (here, k B , A, E c , and T 
are the Boltzmann constant, gap, charging energy, and 
temperature, respectively), only two charge states, n = 
and n = 1, are important for the dynamics of the sys- 
tem, and thus this device behaves as a two-level sys- 
tem {|0) , 1 1)}. For N identical charge qubits placed in a 
single-mode cavity 1|| 2(J (Fig. 1(d)), one can select a 
subset of qubits for the gate, while the remaining qubits 
are not involved in the gate operation, by setting their 
$ = $ /2, U/ c = e/C g , and V g &c = to have them de- 
coupled from the cavity mode and their free Hamiltonian 
being zero. Here, $0 is the flux quantum, <E> is the exter- 
nal magnetic flux piercing the SQUID loop, V g c is the 
dc gate voltage, V g c is the ac gate voltage, and C g is 
the gate capacitance. The method presented below for a 
NTCP gate works for a subset of qubits selected from the 
N qubits in a cavity, because the qubit-qubit coupling, 
mediated by the cavity mode, does not depend on the 
relative position between any two qubits. Without loss 
of generality, we assume that the set of qubits involved 
in the gate operation are the n + 1 qubits labelled by 
1, 2, and n + 1 (here, 1 < n < N). The Hamiltonian 
for the n+1 qubits and the cavity mode is 

H = hu c a^a + E z (V g dc ) S z - Ej(<f>) S x 

+httcos(ojt + if) S z + Kg (a + a)) S z , (1) 

where a is the photon annihilation operator of the cavity 
mode with frequency co c ; S z and S x are the collective op- 
erators of the qubits, given by S z — Y^j=i a z,j an d S x — 



YTj=l a x,j, with Pauli operators a z j = (O^. (0| - (1| 
and a Xt j = \0). (1| + |1) . (0| for qubit j; g is the coupling 
constant between the cavity mode and each qubit; and Cl 
is the Rabi frequency. In addition, E z = —2E C (1 — 2n g c ) 
with the charge energy E c = e 2 /(2C g + 4C,/o) and 
n g c = C g V g c j (2e) . The effective Josephson coupling is 
given by Ej ($) = 2E J0 cos (7r$/$ ) ■ The fourth term 
of Eq. (1) comes from the ac gate voltage given by 
y-ac _ y o cog ^-j- _|_ ^ w hile the last term arises from 
the quantum part of the gate voltage given by U ff qu = 
V^ qu (o + a*), which is caused by the electric field of the 
cavity mode when a qubit is inside the cavity. Finally, 
the coupling constant g and the Rabi frequency £1 are 
given by g = 2E c C g V<? u / {he) and ft = 2E c C g V / (he) , 
respectively. Hereafter, we set E z = (i.e., n g c = 1/2), 
ujq = Ej (<E>) /h, and h = 1. From the last two terms 
of Eq. (1), in the interaction picture with respect to 
Hq = uj c a^a — ojqS x , we obtain (under a rotating-wave 
approximation and assuming lu = 2lvq): 



r 

Hi = — [S z cos ip 



i{S + - S-)smip] , (2) 
a(S z + S- -S+) +H.c] , (3) 



where S + and S are, respectively, the raising and lowing 
operators for the qubits, given by S + — (0| 
and S~ — Y^j—i |0)j (1| i and 5 = lo c — uj is the detuning 
between the cavity-mode frequency and the frequency of 
the ac gate voltage. 

Unitary evolution. — We now consider two special 
cases: <p = and 8 < 0, as well as ip = ir and 5 > 0. The 
results from the unitary evolution, obtained for these two 
special cases, will be employed below for the gate imple- 
mentation. Let us begin with the case (p = and i5 < 0. 
Performing the transformation e~ tHlt on H2 for ip = 0, 
and assuming fl 3> \S\ , g, we obtain 



H' 2 = e lUlt H 2 e 



Hit tt --Wit _ 9 f—iSt 



v 1 —tot 1 iSt +\ o 

- (e a + e a ] )b z 



(4) 



The evolution operator U' for H' 2 takes the form [13j, 



JJ' (t) = e -iMt)Sl e -iB(fiS x a e -iB*(t)S,at /g\ 

with B{t) = ig(e- lSt - l)/(25) and A{t) = g[2B*(t) - 
gt]/(45). When t = r = 2tt/ \5\ , we have B(t) = and 
A(t) = —g 2 r /(4(5). Then, the evolution operator (in the 
Schrodinger picture) of the qubit system is 



U (t) = e 



—IHqT —iH\TTTt 



U'(r) 



= cxp^woTS'z) exp(— iQrSz/2) exp(— i\rS z ), (6) 

where A = — g 2 / (45) > 0. Note that for qjqt = rrni 
(hereafter m is an interger), one has exp(iujoTS x ) — 
(— l) m 11^=1 hi w here Ij is the identity operator for 
qubit j. Thus, the operator U reduces to 

U (r) = exp(-iClTS z /2) exp(-i\rS 2 ). (7) 
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Let us now consider the case ip = ir and 5 > 0. For con- 
venience, we replace 8 and g by 8' and g' , respectively. By 
setting E z — Ej ($) = Q = for qubit 1, the coupling of 
qubit 1 with the cavity mode becomes negligibly small. 
The terms corresponding to j = 1 can thus be dropped off 
from H, Hq , Hi , and Hi . Note that the reduced Hamil- 
tonian ([2]) and ([3]) for tp = tt have forms similar to those 
for tp = 0. Therefore, it is straightforward to show that 
under the condition ^> 8', g', when t = r' = 2tt/8', the 
evolution operator of the qubits becomes 

U (t) = exp{iu T'S' x ) exp(irLT'S' z /2) expiiX'r'Sf), (8) 

where A' = g' 2 / (AS') , S z = v.j and S£ = 

E^=2 °"a!)j- The operator U can be reduced to 

U (r') = exp(«nr'5^/2) exp{iX' t' S' z 2 ) (9) 

for o^ot' = to7t. Finally, we note that the operator U (or 
U) does not include the photon operator a or a + of the 
cavity mode. Hence, the cavity mode can be initially in 
an arbitrary state (e.g., in a vacuum state, a Fock state, 
a coherent state, or even a thermal state). 

Implementation of a NTCP gate. — The operations 
for the gate implementation and the unitary evolutions 
after each step of operation are listed below: 

Step (i): Set V dc = e/C g and V™ = V cos (cot) for 
each qubit. Choose $ (applied to each qubit) and uj ap- 
propriately such that u>o($>) = w/2 = 0.5moj c / (m — 1) , 
leading to 8 = —uj c / (m — 1) < (m > 1) and thus 
uiqt = rn7r for r = —2tt/8. One can see that this is the 
case discussed above for tp — and 8 < 0. Thus, the U 
of Eq. ([7} is the evolution operator for the qubit system 
for the time r = —2tt/S. 

Step (ii): Set V* c = e/C g , V* c = 0, and $ = $ /2 for 
qubit 1; while set V dc = e/C* 9 and V™ = V cos (wt + tt) 
for qubits (2, 3, n + 1). In addition, set uiq(&) — 
lu/2 = 0.5 (m- 2)w c / (m - 1) for qubits (2,3, ...,n+ 1), 
by choosing $ and a; appropriately. Accordingly, we 
have 5' = \S\ — lu c / (to — 1) > (m > 2) and thus 
ujqt' = (m — 2) 7r for r' = 2tt/5'. One can see that this is 
the case discussed above for tp = it and 8 > 0. Hence, for 
an evolution time r' = 2ir/5', the evolution operator of 
the qubit system is f in Eq. (J9j> . 

When (5' = —8 and g' = g, we have A' = A and t' = r. 
In this case, the joint time evolution operator, after the 
above two-step operation, is 

U(2t) = U{t)U{t) 

= cxp (— if2rcr Z) i/2) exp (— i2\Tu Zy iS' z ) . 

The condition = — 5 is automatically satisfied by the 
steps above. 

Step (iii): Set V* c = and $ = $ /2 for each qubit. 
Set V g dc = 2enf }1 /C g for qubit 1 while V dc = 2enf/C g 
for qubits (2,3, ...,n + 1). In addition, adjust the cavity- 
mode frequency [i^, such that it is highly detunned 



with the |0) «-> |1) transition of each qubit. Thus, the 
Hamiltonian in Eq. (1) for the qubit system becomes 
H = E Zil a z<1 + E Z S' Z , where E z>1 = -2E c (l - 2n^) 
while E z — —2E C (1 — 2n dc ). For a time interval r, the 
corresponding evolution operator is then given by 

U (r) = exp(— iE Z! iTa Zi i/h) exp(— iE z rS' z /K). 

With a choice of n dc A = 0.5 - h{An\ + fl)/(8E c ) 
and n dc = 0.5 - h\/(2E c ), one can find from U(2t) 
and U (t) above that the joint time evolution oper- 
ator, after the above three-step operation, is given 
by C/(3r) = U(t)U(2t) = l^^U p (l,j), with 
U p (l,j) = ex.p[i2\r (a z .i + a z j - a z .ia z .j)}. For the 
qubit pair (l,j), we have U p (l, j) |ri) \sj) = |n) \sj) 
(with rs = 00, 01, or 10), while U p (1, j) |li) ] Xj> = 
cxp(— «8At) 1 1 1 > |Lj) , where an overall phase factor 
exp(z2Ar) is omitted. This result shows that for 8Ar = 
(2fc + l)7r, i.e., \S\ — 2gj\j2k + 1 (fc is an integer), a 
two-qubit controllcd-phase gate described by U p (1, j) = 
Iij — 2 1 1 j 1 j ) (lilj| is achieved for the qubit pair (1, j). 
Here and below, qubit 1 acts as a control while qubit j 
as a target, and I\j is the identity operator for the two 
qubits 1 and j. All above conditions on 8 can be satisfied 
with an appropriate choice of m, k, lu c and g. 

Finally, we have U(3r) = U^ihj ~ 2 |lilj) (lil,-]), 
which demonstrates that n two-qubit controlled-phase 
gates are simultaneously performed on the qubit pairs 
(1,2), (1, 3),..., and (1, n+1), respectively. Note that each 
qubit pair contains the same control qubit (i.e., qubit 1) 
and a different target qubit. Hence, a NTCP gate with n 
target qubits (2,3, ...,n + 1) and one control qubit (i.e., 
qubit 1) is obtained after the above three-step operation. 
A brief overview on this NTCP gate is in Fig. 1(a). 

It can be found from U(3t) and U p (l,j) above that 
the present method is based on an effective Hamiltonian 

n+l n+1 

H cS = y^Hij = -2A^ (a z> i + a z j - a Z: ia zJ ) , (10) 

j=2 3=2 

which contains the interaction terms between the con- 
trol qubit and each target qubit, but does not include 
interaction terms between any two target qubits. Note 
that each term Hu in Eq. (JTUJ) acts on a different tar- 
get qubit, with the same controlled qubit, and that any 
two terms H\j for different j commute with each other. 
Therefore, the n two-qubit controlled-phase gates form- 
ing the NTCP gate can be simultaneously performed on 
the qubit pairs (1, 2), (1, 3),..., and (1, n + 1). 

Proposed experiment. — For this method to work, the 
following conditions should be met. For steps (i) and 
(ii), the Rabi frequency f2 needs to be much larger than 
<?, g' , — 8, 8'; the deviation from the degeneracy point 
is given by So — H (fi + g) j (4E C ) , which needs to be a 
small number to have the qubits working near the degen- 
eracy point. For step (iii), from n d \ and n dc above, we 
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FIG. 2: (Color online) Proposed setup for six charge qubits 
(red squares) and a (grey) standing-wave quasi-one dimen- 
sional coplanar waveguide cavity. Each qubit is placed at an 
antinode of the electric field. The two blue curves represent 
the standing-wave electric field, along the y-direction. 

obtain e 1 = h(4n\ + tt)/(8E c ) and e 2 = HX/(2E C ). Here, 
£\ is the deviation from the degeneracy point for the con- 
trol qubit 1 while £i for the n target qubits, which should 
be sufficiently small to ensure that the qubits work near 
their degeneracy points. 

As a concrete example, let us consider the experimental 
feasibility of implementing a five-target-qubit controlled- 
phase gate using superconducting charge qubits with 
C g = 1 aF, Cj Q = 300 aF, E c /h = 32 GHz, E J0 /h = 5 
GHz, T-2 = 500 ns, and T x = 7.3 fjs. The charge 
qubits with these parameters are available at present [25l . 
|26[. For a superconducting ID standing- wave coplanar 
waveguide cavity and each qubit placed at an antin- 
ode of the cavity field (as shown in Fig. 2), the ampli- 
tude of the quantum part of the gate voltage is given 
by [13] Vq U = (hujc) 1 ^ 2 (Lcq)' 1 ^ 2 , where L is the cav- 
ity length and cq is the capacitance per unit length 
of the cavity. The coupling constant is then given by 
g = 2E c C g (he) -1 (hujc) 1 ^ 2 (Lcq) -1 ^ 2 , showing that g 
does not depend on the detuning S. Therefore, the con- 
dition g = g' required above can be satisfied. For charge 
qubits with the above parameters and a cavity with 
uj c / (2tt) = 10 GHz, L = A ~ 12 mm, c ~ 0.22 aF/^m, 
and s e = 6.3, a simple calculation giv es q /27r ~ 100 
MHz, which is experimentally available [26j]. Here, A is 
the wavelength of the cavity mode and e e is the effective 
relative dielectric constant. With the choice of \S\ ~ 0.9g 
(corresponding to the integers m — 112, k — 2), the to- 
tal operation time t op = 3r would be ~ 33 ns, which is 
much shorter than the dephasing time T% and the cavity- 
mode lifetime kT 1 — Q/ui c ~ 159 ns for a cavity with 
Q = 10 4 . Here, Q is the (loaded) quality factor of the 
cavity. Note that a quality factor Q = 10 4 has been 
demonstrated by cavity QED experiments with supercon- 
ducting charge qubits [27J ■ For a qubit-cavity system with 
the parameters given above, we have e$ ~ 5.46 x 10~ 3 , 
Ei ~ 4.51 x 10~ 3 , and e 2 ~ 4.34 x 10~ 4 for fi/ (2tt) - 600 
MHz. Therefore, the conditions for the qubits to work 
near the degeneracy point are well satisfied. 



Note also that the adjustment of the cavity-mode fre- 
quency is unnecessary because one can adjust the dc gate 
voltage to have the qubits decoupled from the cavity 
mode. However, when a large number of qubits are in- 
volved, it is much more convenient to adjust the cavity 
mode frequency rather than adjusting the dc gate voltage 
of each qubit. This is because all qubits can be simulta- 
neously decoupled from the cavity mode by adjusting the 
cavity mode frequency, but one will need to individually 
adjust step-by-step the dc gate voltage for each qubit to 
decouple the qubits from the cavity mode. 
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